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PREFACE 


The  study  described  herein  was  performed  at  the  US  Army  Engineer  Water¬ 
ways  Experiment  Station  (WES)  during  February-April  1991  for  the  Headquarters, 
US  Army  Corps  of  Engineers  (HQUSACE) ,  as  part  of  the  Civil  Works  Research  and 
Development  Program.  Funds  were  allocated  under  Civil  Works  Investigation 
Work  Unit  32541,  "Riprap  Design  and  Cost  Reduction;  Studies  in  Near  Prototype 
Size  Laboratory  Channel,"  under  HQUSACE  Program  Monitor  Mr.  Thomas  Munsey. 

This  study  was  accomplished  under  the  direction  of  Messrs.  F.  A. 

Herrmann,  Jr.,  Director  of  the  Hydraulics  Laboratory  (HL) ;  R.  A.  Sager, 
Assistant  Director,  HL;  and  G.  A.  Pickering,  Chief  of  the  Hydraulic  Structures 
Division  (HSD) ,  HL.  The  tests  were  conducted  by  Dr.  S.  T.  Maynord,  project 
engineer,  and  Mr.  D.  M.  White,  Spillways  and  Channels  Branch  (SCB),  HSD,  under 
the  direct  supervision  of  Mr.  N.  R.  Oswalt,  Chief,  SCB.  This  report  was 
written  by  Dr.  Ma3mord  and  edited  by  Mrs.  M.  C.  Gay,  Information  Technology 
Laboratory ,  WES . 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Leonard  G.  Hassell,  EN. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used 
(metric)  units  as  follows: 

Multiply 
cubic  feet 
degrees  (angle) 
feet 
inches 


in  this  report  can  be  converted 

_ Ex _ 

0.02831685 

0.01745329 

0.3048 

2.54 


to  SI 

To  Obtain 
cubic  metres 
radians 
metres 
centimetres 
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RIPRAP  RESISTANCE  TESTS  FROM  A  LARGE  TEST  CHANNEL 


PART  I :  INTRODUCTION 

Background 

1.  Resistance  of  flow  over  riprap  boundaries  Is  an  Important  part  of 
channel  design  and  Is  used  In  determining  water— surface  elevations  and  veloci¬ 
ties.  A  previous  study*  dealing  with  the  flow  resistance  of  riprap  used  data 
from  several  rectangular  tilting  flumes  In  which  the  riprap  was  placed  on  the 
bottom  of  the  flume.  Results  from  that  study  showed  that  flow  resistance  for 
relative  depth  d/Ogg  with  values  from  3  to  30  Is  best  described  by  a  power 
law  rather  than  the  commonly  used  logarithmic  function.  Here  d  Is  depth  and 
Dgo  Is  the  particle  size  of  which  90  percent  Is  finer  by  weight.  Strlckler's 
equation**  was  derived  from  the  power  law  equation  and  resulted  In  the 
following,  using  particle  size  Dgg  In  feet 

n  *  0.0360  Dgo^'®  d) 


and  using  particle  size  D^g  In  feet 


n  =  0.0380  Djgi'® 


(2) 


where 

n 

Dso 


Manning  roughness  coefficient 

particle  size  of  which  50  percent  Is  finer  by  weight 


*  Steve  T.  Maynord.  1991.  "Flow  Resistance  of  Riprap,"  Journal  of 
Hydraulic  Engineering.  American  Society  of  Civil  Engineers.  Vol  117, 

No.  6,  pp  687-696. 

**  A.  Strickler.  1923.  "Contributions  to  the  Question  Concerning  a  Formula 
for  Speed  and  the  Roughness  Numbers  for  Rivers,  Channels  and  Culverts" 
("Beltrage  zur  Frage  der  Geschwindigkeicformel  und  der  Rauhigkeitszahlen 
fur  Strome,  Kanule  und  Geschlossene  Leltungen"),  Mlttellungen  des  Amtes  fur 
Wasserwlrtschaft,  No.  16,  Bern,  Switzerland,  pp  12-13  (in  German). 
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Equations  1  and  2  are  based  on  data  from  wide  rectangular  flumes  having 
essentially  two-dimensional  flow. 


Purpose  and  Scope 

2.  The  objective  of  this  study  is  to  develop  techniques  for  estimating 
flow  resistance  over  riprap  boundaries.  The  scope  of  this  report  is  to  com¬ 
pare  flow  resistance  data  collected  in  the  trapezoidal  channel  in  the  Riprap 
Test  Facility  (RTF)  at  the  US  Army  Engineer  Waterways  Experiment  Station  to 
the  results  of  the  data  analysis  from  the  rectangular  tilting  flume  study. 
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PART  II:  EXPERIMENTAL  INVESTIGATION 


3.  The  RTF  (Plate  1)  is  a  large  outdoor  test  channel  having  a  780— ft* 
length,  four  bendways,  and  a  discharge  range  of  0  to  200  cfs.  Two  constant- 
speed  pumps,  Cl  and  C2,  provide  discharges  of  52  and  49  cfs,  respectively. 

Two  variable-speed  pumps ,  VI  and  V2 ,  provide  maximian  discharges  of  42.5  and 
48  cfs,  respectively.  The  initial  203  ft  of  the  RTF  (Plate  2)  was  used  in 
this  study,  which  is  a  straight  reach  having  a  trapezoidal  section  with  12-ft 
bottom  width.  The  channel  has  1V:2H  side  slopes  on  both  sides  to  sta  1+71 
followed  by  an  8-ft-long  transition  on  the  left  descending  bank  to  a  side 
slope  of  1V:1.5H.  The  first  20  ft  of  the  channel  is  covered  with  6-  to 
12-in. -diam  riprap  to  dissipate  turbulence  of  flow  leaving  the  pump  discharge 
flume.  The  remainder  of  the  straight  reach  was  covered  with  riprap  having  the 
gradation  shown  in  Plate  3.  Plastic  tubing  was  placed  beneath  the  riprap  at 
sta  0+47,  0+76,  1+05,  1+34,  and  1+63  and  was  connected  to  a  stilling  well. 

The  sensing  end  of  the  tube  terminated  at  the  channel  center  line.  The 
channel  cross  section  was  surveyed  at  11  sections  and  results  are  shown  in 
Table  1.  The  elevations  of  the  five  survey  points  on  the  channel  bottom  were 
averaged  at  each  cross  section  and  are  plotted  in  Plate  4.  A  least  squares 
fit  of  the  average  bottom  elevation  data  resulted  in  a  bottom  slope  of 
0.00289  ft/ft  from  sta  0+47  to  1+71. 

4.  Vertical  velocity  profiles  were  determined  at  sta  0+47,  1+05,  and 
1+63  to  determine  if  the  velocity  profile  was  fully  developed  at  sta  0+47, 
which  was  27  ft  from  the  large  change  in  boundary  roughness  at  sta  0+20.  If 
the  profile  was  not  fully  developed,  data  from  sta  0+47  might  be  in  error. 
Profiles  were  measured  at  the  channel  center  line  and  2  ft  from  each  side  of 
the  center  line.  Results  are  shown  in  Plates  5-7  and  8-10  for  discharges  of 
52  and  76  cfs,  respectively.  The  velocity  profile  at  sta  0+47  is  not  signifi¬ 
cantly  different  from  those  of  the  downstream  stations,  which  means  that  data 
from  sta  0+47  can  be  used  in  the  analysis. 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  metric 
(SI)  units  is  found  on  page  3. 
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PART  III:  ANALYSIS  AND  RESULTS 


5 .  Water-surface  elevations  were  measured  for  nine  tests  as  shown  in 
Table  2  and  plotted  in  Plates  11  and  12.  Note  that  the  75—  and  76-cfs  pro¬ 
files  exhibit  a  marked  drop  in  water  surface  between  sta  1+34  and  1+63.  This 
is  likely  an  effect  of  the  downstream  transition  that  was  not  significant  at 
the  lower  discharges.  Comparisons  of  observed  versus  computed  water-surface 
elevations  will  not  use  sta  1+63  data  for  75  and  76  cfs.  All  other  profiles 
are  typical  backwater  curves.  Tests  1-8  represent  dry-placed  riprap  having  a 
relatively  even  surface.  Test  9  ^  as  conducted  with  the  riprap  roughened  to 
simulate  the  rougher  surface  t3rpical  of  riprap  placed  underwater. 

6.  The  cross-section  data  in  Table  1  were  used  in  the  HEC-2  water — 
surface  profile  computation.  For  discharges  other  than  75  and  76  cfs,  the 
observed  water-surface  elevation  at  sta  1+63  was  used  as  the  downstream  start¬ 
ing  water-surface  elevation  in  the  computations.  For  discharges  of  75  and 

76  cfs,  the  observed  water-surface  elevation  at  sta  1+34  was  used  as  the  down¬ 
stream  starting  water-surface  elevation.  The  n  value  that  resulted  in  the 
best  agreement  between  observed  and  computed  water-surface  elevation  is  shown 
in  Table  2  along  with  the  corresponding  computed  water-surface  elevations , 

The  average  n  value  for  Tests  1-8  was  0.0272.  Using  D90  “  1.93  in./12  - 
0.161  ft  from  Plate  3  and  the  Strickler  equation*  results  in 

n  -  0.0369  (3) 


Using  D50  -  1.4  in./12  -  0.117  ft  results  in 


n  =  0.0389 


(4) 


The  coefficients  in  Equations  3  and  4  are  about  2.4  percent  greater  than  Equa¬ 
tions  1  and  2  from  Maynord. 

7.  The  logarithmic  equation  for  flow  resistance  is 


*  Strickler,  op.  cit. 
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(5) 


5.75  log 


where 

a  -  channel  shape  factor  given  by  Hey* 

R  -  hydraulic  radius 
Kg  -  equivalent  sand  grain  roughness 

and  f  is  the  Darcy  friction  factor  defined  as 


f  = 


(6) 


where 

g  -  acceleration  due  to  gravity 
S  -  channel  slope 
V  -  channel  velocity 

The  previous  study  found  that  the  logarithmic  equation  (5)  is  not  valid  for 
intermediate  scale  roughness  having  3  ^  R/Dgo  ^  10.  The  data  in  Table  3  were 
used  to  determine  the  least  squares  logarithmic  equation 

(ir-'‘-'0  1og|^j*6.29  (7) 


As  in  the  previous  study,  Student's  t-test  shows  that  the  slope  of  Equation  7 
of  4.20  is  statistically  different  from  the  commonly  accepted  value  of  5.75. 
Because  of  this  difference,  this  study  suggests  that  the  logarithmic  equations 
are  not  valid  over  the  range  of  data  used  herein  having  4.5  <  R/Dgo  <  8.2. 

8.  Comparison  of  results  from  Test  9  with  the  roughened  riprap  to 
Test  2  with  the  dry-placed  riprap  (Table  2)  shows  a  (0. 0306-0. 0270)/0. 0270 
-  13  percent  increase  in  n  for  the  roughened  riprap  that  was  intended  to 
simulate  riprap  placed  underwater. 


*  R.  D.  Hey.  1979.  "Flow  Resistance  in  Gravel-bed  Rivers,"  Journal  of  the 
Hydraulics  Division.  American  Society  of  Civil  Engineers.  Vol  105,  No.  HY4, 
pp  365-379. 
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PART  IV;  CONCLUSIONS  AND  RECOMMENDATIONS 


9.  Results  from  resistance  tests  in  the  RTF  are  in  agreement  with 
results  from  tests  in  the  rectangular  tilting  flvunes.*  The  coefficient  in 
Strickler's  equation**  from  the  RTF  was  2.4  percent  greater  than  the  coeffi¬ 
cient  for  the  rectangular  tilting  flume.*  Results  from  this  study,  similar  to 
the  previous  study,*  suggest  that  the  logarithmic  equation  is  not  valid  for 
intermediate  scale  roughness . 

10.  Riprap  placed  to  simulate  underwater  placement  had  an  n  value 
13  percent  greater  than  riprap  placed  in  dry  conditions,  where  the  rock 
surface  is  much  smoother. 


*  Maynord,  op.  cit. 

**  Strickler,  op.  cit. 
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Table  1 

Measured  Cross  Sections 


Station  0-f30  Station  0+47  _ Station  0+61,^ 


acai 

X 

LJ.U11  VTJV 

Elevation 

X 

Elevation 

X 

Elevation 

6.40 

101.91 

6.41 

101.98 

6.35 

101.83 

8.22 

101.05 

8.25 

101.17 

8.18 

101.01 

10.00 

100.19 

10.00 

100.21 

10.00 

100.19 

13.00 

100 . 09 

13.05 

100.10 

13.00 

100 . 04 

16.00 

100.09 

16.10 

100.11 

16.00 

100.08 

19.00 

100.15 

19.15 

100.07 

19.00 

100.14 

21.75 

100 . 20 

22.20 

100.24 

21.80 

100.22 

23.60 

100.95 

23.98 

101.16 

23.62 

101.06 

25.39 

101.85 

25.83 

101.92 

25.47 

101.81 

Station  0+76 

Station  0+90.5 

Station  1+05 

X 

Elevation 

X 

Elevation 

X 

Elevation 

6.39 

101.87 

6.38 

101.79 

6.38 

101.78 

8.19 

101.00 

8.19 

100.94 

8.20 

100.95 

10.00 

100 . 14 

10.00 

100.09 

10.00 

100.08 

13.00 

100.03 

13.00 

99.94 

13.00 

100.00 

16.00 

99.99 

16.00 

99.96 

16.00 

99.93 

19.00 

100.02 

19.00 

100.00 

19.00 

100.00 

22.00 

100.07 

21.80 

99.99 

22.00 

100.08 

23.82 

100.90 

23.59 

100.80 

23.80 

100.95 

25.59 

101.83 

25.39 

101.75 

25.63 

101.76 

Station  1+19.5 

Station  1+34 

Station  1+48 . 5 

X 

Elevation 

X 

Elevation 

X 

Elevation 

6.39 

101,77 

6.43 

101.75 

6.35 

101.59 

8.18 

100.88 

8.23 

100.87 

8.13 

100.68 

10,00 

100 . 04 

10.00 

99.94 

10.00 

99.96 

13.00 

99.93 

13.00 

99.83 

13.00 

99.80 

16.00 

99.88 

16.00 

99.80 

16.00 

99.79 

19.00 

99.87 

19.00 

99.80 

19.00 

99.82 

22.00 

99.97 

21.80 

99.91 

21.75 

99.87 

23.79 

100.87 

23.65 

100.68 

23.56 

100.73 

25.59 

101.74 

25.48 

101.47 

25.32 

101.68 

Station 

1+63 

Station  1+71 

X  Elevation 

X 

Elevation 

6.39 

101.64 

6.41 

101.68 

8.16 

100.70 

8.19 

100.78 

10.00 

99.92 

10.00 

99.92 

13.00 

99.79 

13.00 

99.75 

16.00 

99.77 

16.00 

99.75 

19.00 

99.80 

19.00 

99.74 

22.00 

99.89 

22.00 

99.82 

23.79 

100.78 

23.80 

100.69 

25.63 

101.58 

25.61 

101.55 

Note : 

Elevations  are 

referenced  to  an 

arbitrary 

datum. 

X  -  distance  from  left  descending  side  of  channel,  feet. 


Table  2 

Observed  Versus  Computed  Water-Surface  Elevations 


Test 

Discharge 

cfs 

Best 

Uater-Surface  Elevation*  at  Sta 

f«o, 

n 

Pumns 

0+47 

0+76 

1+05 

1+34 

1+63 

1 

52.0 

0.0278 

Cl 

101.480 

101.480 

101.421 

101.420 

101.368 

101.360 

101.326 

101.330 

101.281 

101.280 

2 

48.0 

0.0270 

V2 

101.394 

101.390 

101.342 

101.340 

101.286 

101.280 

101.246 

101.250 

101.199 

101.200 

3 

49.0 

0.0272 

C2 

101.423 

101.420 

101.361 

101.370 

101.310 

101.300 

101.271 

101.270 

101.227 

101.230 

4 

42.5 

0.0268 

VI 

101.299 

101.300 

101.236 

101.240 

101.183 

101.180 

101.143 

101.150 

101.098 

101.100 

5 

25.0 

0.0277 

V2 

100.998 

101.000 

100.929 

100.930 

100.872 

100.870 

100.830 

100.830 

100.788 

100.790 

6 

25.0 

0.0273 

VI 

100.987 

100.990 

100.918 

100.920 

100.861 

100.860 

100.818 

100.820 

100.776 

100.780 

7 

76.0 

0.0268 

Cl,  V2 

101.797 

101.790 

101 . 740 
101.750 

101.688 

101.690 

101.660 

101.660 

101.566 

8 

75.0 

0.0272 

C2,  VI 

101.783 

101.780 

101.727 

101.730 

101.678 

101.670 

101.643 

101.640 

101.548 

9 

48.0 

Disturbed 

0.0306 

V2 

101.441 

101.430 

101.358 

101.370 

101.295 

101.290 

101.241 

101.240 

101.175 

101.180 

*  Elevations  are  referenced  to  an  arbitrary  datum.  The  first  row  of  data  fot 
each  test  gives  the  observed  elevations,  and  the  second  row  gives  the 
computed  elevations 


Hydraulic 


Test 

Depth* 

Radius 

No. 

ft 

R.  ft 

1 

1.35 

1.10 

2 

1.27 

1.04 

3 

1.29 

1.06 

4 

1.16 

0.97 

5 

0.85 

0.74 

6 

0.84 

0.73 

7 

1.67 

1.31 

8 

1.66 

1.31 

9 

1.28 

1.05 

n 


_S_ 

(from 

fs! 

D90 

Table  2) 

UJ 

6.84 

0.0278 

9.57 

6.47 

0.0270 

9.76 

6.59 

0.0272 

9.72 

6.03 

0.0268 

9.72 

4.60 

0.0277 

8.99 

4.54 

0.0273 

9.10 

8.15 

0.0268 

10.22 

8.15 

0.0272 

10.07 

6.53 

0.0306 

...I 

*  Water-surface  elevation  -  average  bottom  elevation  at  sta  1+05. 


RIPRAP  TEST  FACIUTY  SCHEMATIC 


PLATE  3 


SIEVE  SIZE 


BOTTOM  ELEVATION.  FT 


Y/D 


LEGEND- 

Y  =  DISTANCE  ABOVE  BOTTOM 
D  =  DEPTH 

V  =  VELOCITY  AT  Y 

V^VG  =  DEPTH-AVERAGED  VELOCITY 


NOTE- 

VELOCITY  PROFILES  AT 
CHANNEL  CENTER  LINE  AND 
2  FT  ON  EACH  SIDE,  CURVES 
NOT  LABELED  DUE  TO  SIMILARITY 


DIMENSIONLESS  VELOCTTY 
PRORLES 

DISCHARGE  52  CFS 
STATION  1+05 


PLATE  6 


Y/D 


v/Vavg 


LEGENDi 

Y  =  DISTANCE  ABOVE  BOTTOM 
D  =  DEPTH 

V  =  VELOCITY  AT  Y 

vavg  =  depth-averaged  velocity 


NOTE.  DIMENSIONLESS  VELOCTTY 

VELOCITY  PROFILES  AT  PROFILES 

CHANNEL  CENTER  LINE  AND 

2  FT  ON  EACH  SIDE.  CURVES  DISCHARGE  76  CFS 

NOT  LABELED  DUE  TO  SIMILARITY  STATION  0+47 


PLATES 


Y/D 


v/Vavg 


LEGEND: 

Y  =  DISTANCE  ABOVE  BOTTOM 
D  =  DEPTH 

V  =  VELOCITY  AT  Y 

v^vG  =  depth-averaged  velocity 


NOTE:  DIMENSIONLESS  VELOCTTY 

VELOCITY  PROFILES  AT  PRORLES 

CHANNEL  CENTER  LINE  AND 

2  FT  ON  EACH  SIDE,  CURVES  DISCHARGE  52  CPS 

NOT  LABELED  DUE  TO  SIMILARITY  STATION  1+63 


PLATE? 


Y/D 


0.00  0.25  0.50  0.75  1.00  1.25  1.50 


v/Vavg 

LEGEND! 

Y  =  DISTANCE  ABOVE  BOTTOM 
D  =  DEPTH 

V  =  VELOCITY  AT  Y 

VaVG  =  DEPTH-AVERAGED  VELOCITY 


NOTE.  DIMENSIONLESS  VELOCfTY 

VELOCITY  PROFILES  AT  PRORLES 

CHANNEL  CENTER  LINE  AND 

2  FT  ON  EACH  SIDE.  CURVES  DISCHARGE  76  CFS 

NOT  LABELED  DUE  TO  SIMILARITY  STATION  1+63 


PUTE  10 


Y/D 


V/VftVG 


LEGEND' 

Y  =  DISTANCE  ABOVE  BOTTOM 
D  =  DEPTH 

V  =  VELOCITY  AT  Y 

V/:^VG  =  DEPTH-AVERAGED  VELOCITY 


NOTE.  DIMENSIONLESS  VELOCTTY 

VELOCITY  PROFILES  AT  PRORLES 

CHANNEL  CENTER  LINE  AND 

2  FT  ON  EACH  SIDE.  CURVES  DISCHARGE  76  CFS 

NOT  LABELED  DUE  TO  SIMILARITY  STATION  1+05 


PLATES 


WATER-SURFACE  ELEVATION,  FT 


102.00 


STATION,  FT 


NOTE-  25-CrS  WATER-SURFACE  PROFILE 

TAKEN  WITH  PUMP  V2 

ELEVATIONS  ARE  REFERENCED  TO 
AN  ARBITRARY  DATUM 


WATER-SURFACE  PROFILES 

DISCHARGE  25,  48,  52,  AND  76  CFS 


PLATE  11 


WATER-SURFACE  ELEVATION,  FT 


NDTEi  S5-CFS  WATER-SURFACE  PROFILE 
TAKEN  WITH  PUMP  VI, 

48-CFS  WATER-SURFACE  PROFILE 
TAKEN  WITH  RIPRAP  SURFACE  ROUGHENED 
TO  SIMULATE  UNDERWATER  PLACEMENT. 
ELEVATIONS  ARE  REFERENCED  TO 
AN  ARBITRARY  DATUM 


WATER-SURFACE  PROFILES 

DISCHARGE  25,  42,5,  48,  49,  AND  75  CFS 


PUTE  12 
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